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Abstract. We prove that the volumes determined by the lengths of the non-zero 
vectors ±x in a random lattice L of covolume 1 define a stochastic process that, 
as the dimension n tends to infinity, converges weakly to a Poisson process on the 
positive real line with intensity i. This generalizes earlier results by Rogers ( [111 
Thm. 3]) and Schmidt (13, Satz 10]). 

1. Introduction 

Let Gn = SL(n, M) and r„ = SL(n, Z). We will be interested in the space Xn = 
^n\Gn considered as the space of lattices of covolume 1. We let /i„ denote the Haar 
measure on Gn, normalized so that it represents the unique right G^i— invariant 
probability measure on the homogeneous space X„. 

For a random lattice L € X„ we study the lengths of the pairs ±x of non-zero 
vectors. Given a lattice L € Xn, we order the non-zero vector lengths in L as 
< £i < ^2 ^ ^3 ^ • • •) where we count the common length of the vectors x and —x 
only once. For j > 1, we define 



V 



_ 



so that Vj is the volume of an n-dimensional ball of radius ij. Finally, for i > 0, we 
let 

Nt{L) := #{j : V, < t}. 

The first result concerning the statistics of the sequence {ij}'^^ is due to Rogers. 
In [TI] he shows that, for fixed t, Nt{L) has a distribution which converges weakly to 
the Poisson distribution with mean | as n ^ oo. Clearly Rogers' result determines 
the limit distribution of each fixed as n — ?• cx). In particular the volume Vi has 
an exponential distribution with expectation value 2 as n — t- oo. In [13] Schmidt 
presents a more direct study of the distribution of li using another method resulting 
in better error bounds as n — t- oo. 

There is a close connection between these results and lattice packings of spheres: 
For each fixed n the random variable Vi has bounded range, i.e. supj;^^^^ Vi < oo. 
The quantity T>n '■= 2""" sup^^gjj^^ Vi is in fact the maximal density of a lattice sphere 
packing in M". The search for upper bounds for the density of lattice (and also more 
general) sphere packings is a very well studied subject, yet optimal bounds are known 
only in very special cases. At present the best known bounds as n ^ oo are due to 
Ball ([!]) and Kabatyanskii and Levenshtein ([6]); these state that 

2(n - l)C(n)2-" < P„ < 2-"(°-599+o(i)) ^ 
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Note that the hmit distribution result for Vi quoted above immediately implies that 
^ oo as n — )• oo. This observation can be improved using knowledge on the 
rate of convergence to the limit distribution; in fact Schmidt can give a lower bound 
for Dn which differs from that of Ball only by a modest constant factor (cf. [1'6\ Thm. 
11]). 

The results of Rogers and Schmidt suggest very naturally that it should be possible 
to understand not only the distribution of each fixed ij but also the joint distribution 
of the entire sequence for a random lattice L G X„ as n ^ oo. The purpose 

of the present note is to point out that this is possible using Rogers' methods. Our 
main theorem states that, as n ^ oo, the volumes {Vj}j%i behave like the points of 
a Poisson process on the positive real line with intensity ^. 

Theorem 1. For any fixed N, the N -dimensional random variable (Vi,...,VAr) 
converges in distribution to the distribution of the first N points of a Poisson process 
on the positive real line with intensity ^ as n —)• oo. 

Or, put equivalently: 

Theorem 1'. Let {N{t),t > 0} be a Poisson process on the positive real line 
with intensity ^. Then the stochastic process {Nt{-),t > 0} converges weakly to 
{N{t),t> 0} asn^ oo. 

The author got interested in this problem during recent investigations of the value 
distribution of the Epstein zeta function. We recall that for s > ^ and L £ Xn the 
Epstein zeta function is defined by 

(1) EiL,s)= ^'|mr'^ 

where ' denotes that the zero vector should be omitted. E{L, s) has an analytic 
continuation to C except for a simple pole at s = ^. Now Theorem [1] together with 
equation ([1]) suggest that for a fixed c > ^ the distribution of the random variable 
E{-,cn) (suitably normalized) converges weakly to the distribution of a similar sum, 
defined in terms of the points of a Poisson process, as n — > oo. This and related 
results will be presented in detail in |14j . 

Let us also point out the close conceptual relation between our Theorem [1] and 
the main result (Theorem 1.1) by VanderKam in [15] (cf. also Sarnak [l2]), which 
states that for a generic fixed lattice L £ Xn, the local spacing statistics of the 
infinite sequence {Vj}'^i exhibit Poissonian behavior (to a larger extent the larger 
the dimension). In precise terms: 

Theorem 2 (VanderKam). For almost all L € Xn the 2-, 3-,..., and [n/2]-level 
correlations of the sequence {Vj}^^ are Poissonian, i.e. for any given intervals 
Ii, . . . , Im-i C M (2 < m < |) we have 

^.^ #{{il,...,in^)e{Z+)T■.V^„...,V^^<N,Vi,-V,^GI,-^} ^ . 

i=i 

where (Z"*")^ is the set of all m-tuples of distinct positive integers. 

VanderKam's theorem is of particular interest in the context of the Berry- Tabor 
conjecture (cf. [151 Cor. 1.2]). This conjecture states that for a "generic" completely 
integrable system the distribution of the local spacings between the eigenvalues of 
the Laplacian (suitably normalized) should match that of the points of a certain 
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Poisson process on the positive real line. VanderKam's theorem gives evidence in 
support of this conjecture, since it implies that for almost all flat tori the distribution 
of the eigenvalues of the Laplacian agree, at least up to a certain number of correla- 
tions, with the distribution of the points of a Poisson process. In this context, let us 
note that our Theorem [T] can be viewed as a kind of "low-eigenvalue- high-dimension" 
analog of the Berry- Tabor conjecture for flat tori. Indeed, Theorem [T] can be refor- 
mulated as follows. (Note that the eigenvalues for the flat torus R"/A are 47r^|^p 
with i belonging to the dual lattice A*; "desymmetrizing" and renormalizing these 
to have mean spacing one we get the sequence {|Vj}^^ for A*.) 

Theorem 1". For a random flat torus M'^/A with A G X„ the distribution of the non- 
zero eigenvalues (normalized to have mean-spacing one) converges to the distribution 
of the points of a Poisson process on the positive real line with intensity 1 as n ^ oo. 

2. Rogers' integration formula 

In this section we fix some notation concerning an integration formula that will 
be the major technical tool in the proof of our main result. 

Let 1 < k < n — 1 and let p : (W^)^ — >■ M be a non-negative Borel measurable 
function. In ^ Rogers considers the integral 



/ ^ p{mi,...,mk)dpniL), 



and shows that it equals a certain (positive) infinite linear combination of integrals 
of p over various linear subspaces of (R")*^. In this note we will be interested in the 
similar integral 

(2) / ^ p{mi,...,mk)dpn{L). 

mi,...,mkeL\{0} 

It follows from Rogers' formula in [9j that the integral in ([2]) equals 
(3) 

/ • • • / p{xi, ...,Xk)dxi... dxk 

+ wvf----— ■■■ ofv— ^- s^—x]dx, dx 

Here the outer sum is over all divisions {v, p) = (i^i, . . . , z^^; pi, . . . , pk-m) of the 
numbers 1, . . . ,k into two sequences i^i, . . . and pi, . . . , pk-m with 1 < m < k — 1, 
satisfying 

1 < Ul < 1^2 < . . . < J^m < k, 
(4) 1 < ^1 < ^2 < • • • < ^fc-m < k, 

i^i ^ Pj, if 1 < i < m, 1 < j < k — m. 

The inner sum in ^ is over all m x A; matrices D, with no column vanishing, with 
integer elements having greatest common divisor equal to 1, and with 

diuj = qSij, i = I,. . . ,m, j = I,. . . ,m, 

di^. =0, \i pj <Vi, i = 1, . . . , m, j = I, . . . ,k - m. 
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We call these matrices (z^, ;u)-admissible. Finally Cj = {ei,q), i = l,...,m, where 
El, . . . ,£m are the elementary divisors of the matrix D. 

Remark 1. It follows from the conditions on the matrices D above and [U Thm. 
14.5.1] that in all cases we have ei = 1. In particular it follows that we always have 

\q q J 
3. The joint moments of Nj{L) 

We will now consider the following situation: Let A; > 1 and fix < Vi < V2 < 
• • • ^ Vfc. Given n > k + 1 and 1 < j < we let Nj{L) denote the number of 
non-zero lattice points of L in the n-ball of volume Vj centered at the origin. The 
goal of this section is to establish the following result concerning the joint moments 
of the functions Nj{L): 

Theorem 3. Let k > 1 and fix < Vi < V2 < ■ ■ ■ < V^. For each division {i',iJ,) 
in Rogers' formula let My denote the number of {v, admissible matrices D, with 
elements dij € {0,±1}, having exactly one non-zero entry in each column. Then 

k k m 

(5) ¥.{J{N,{L)) ^\{V, + Y,Mu,,J{V,. 

j = l j = l i=l 

as n 00, where the sum is over all divisions possible in Rogers' formula. 
Remark 2. It is a straightforward calculation to verify that 



m—l 



i=2 

Hence it is possible to rewrite the result in Theorem [3] in a way not involving any 
reference to matrices. With notation as above we have 

k k m—l m 

j=l m=l l=ui<...<Um<k j=2 i=l 

as n — 7- 00. 

3.1. The proof of Theorem [3l We let Pj, 1 1^ j < k, be the characteristic function 
of the n-ball of volume Vj centered at the origin. Then we can write 



iV,(L)= ^ /5,(m), l<j<k. 



Applying Rogers' formula yields 

k k 



(6) 



X / •■■ / W Pj{^^—Xi]dXi...dXm 



It will turn out that the main contribution comes from the terms where q=\ and the 
matrix D has entries dij G {0, ±1} and exactly one non-zero entry in each column. 
We begin by determining the contribution to ([6]) from these terms. 
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Proposition 1. If I < k < n — 1, then 

TT / pj{x)dx +S2y2 ■■■/ Y[ Pji'^'^v^i) dxi.-.dxm 

k m 

i=l (u,y) j=l 

where the sum over D is taken over all (i^, n)- admissible matrices having entries 
dij S {0, ±1} and exactly one non-zero entry in each column. 

Proof. The matrices D on this form satisfies by definition dj^^. = 1, 1 < J < 
If we furthermore let be such that dx^^^^^ = ±1, 1 < £ < k — m, we find that 

k 



n ^ r\j lit 

/ "'/ TT /^i ( dijXj] dxi... dxm 

J — J- * — 

„ „ m k—m 

= ■■■ I Wp^j i^i) n Pi-n {x\i) dxi... dxr, 



3=1 7=1 

Since the volumes Vj are increasing this integral equals 

„ „ m m 

/ " ' / TiP'^J (^j) dXi... dXm = TT Vuj , 

which gives the desired result. □ 
We next give a bound on the contribution from most of the remaining terms in 



Proposition 2. Let 

d 



Jm- JRn fj-^ V 



Xi I dxi . . . dx 



r- 

and let M{D) be the largest value taken by any determinant of an m x m-minor of 
D. Assume that n > max (m(k — m) + l). Then 

l<m<k-l^ ' 

oo 

{v,i,)q=l D ^ ^ {y^f,) D 

q=l 
M{D)=1 

where the remainder term satisfies 

< R{k) < 2"". 

The implied constant depends on k and the volumes Vi,V2, ■ ■ ■ ,Vk but not on n. 

Proof. This is a trivial adaptation of [lOj, Sec. 9] to our situation. (We bound 
yOi,P2, ••• ,/Ofc-i from above by pfc.) □ 

It is immediate to verify that matrices D with g = 1 and M{D) = 1 have all entries 
dij S {0, ±1}. In particular all the matrices in Proposition [1] are on this form. In the 
remaining part of this section we will discuss the contribution to ([6]) coming from 
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matrices D with q = 1 and M{D) = 1 and at least one column containing more than 
one non-zero entry. Our starting point is the following lemma. 

Lemma 1. For any l<i<j<i<k and y G M" we have 

/ / Pi{xi)pj{X2)p£{±Xi±X2 + y)dXidX2 <2{^)^Vi^. 

Proof. Since 

/ / pi{xi)pj{x2)pe{±xi±X2 + y)dxidx2 

jR" JR" 

< / pi{xi)pi{x2)pe{±xi±X2 + y)dxidx2, 

jR" JR" 

the lemma follows from [11] Lemma 5]. □ 

Following Rogers ( |1H Lemma 6] ) we obtain the final estimate needed in the proof 
of Theorem [31 

Lemma 2. Let k > 3 and fix m satisfying 2 < m < k — 1. Let (z^, p) be a division of 
the numbers 1, . . . ,k satisfying Q with our m. Let D be a {v, p)- admissible matrix 
with q = 1, M{D) = 1 and at least one column containing more than one non-zero 
entry. Let l<£<k — mbe such that the leftmost column with more than one 
non-zero entry is pe and let Ai and A2 be minimal with the property that Ai < A2 
and (iAi,/xf / / dx^,^!,. Then 

n n k in 

(7) / •••/ llp,(Y,d,,x^dx^...dx^m<2{l)'^V^^llV,^, 

where Af = {1, 2, . . . , m} \ {Ai, A2}. 

Proof. With the notation introduced above we get 

/ ■■■/ TT /"j ( dijXj] dxi... dxm 
Jr" Jr" ^ -^-^ ' 

JR" JR" j^j^j \ JR" JR" 

m \ 

( ± a^Ai ± a^Aa + ^ dif,^x^ dxx-^dxx^ j Y{ dxj 



< 

m 

X 

i=A2+l -' jeM 



where we have applied Lemma [T] to the inner integral. □ 

We note that for fixed Vi, . . . ,Vk the right hand side of ([7]) is exponentially smaller 
(in n) than the right hand side of ([5]). We stress that also here the implied constant 
depends on Vi , . . . , 14 . Since there are only finitely many matrices on the form in 
Lemma [2] this, together with Propositions [1] and [21 proves Theorem [3l 
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4. The proof of Theorem [T] 

Let us consider a Poisson process {N{t),t > 0} on the positive real line with 
constant intensity ^. Thus N{t) denotes the number of points falling in the interval 
(0, i]. We recall that N{t) is Poisson distributed with expectation value |. We 
further let Ti , T2 , Ts , . . . denote the points of the process ordered in such a way that 
< Ti < r2 < Ts < ■ ■ ■ . 

Proposition 3. Let k > 1 and denote by V{k) the set of partitions of {1, . . . ,A;}. 
For 1 < j < k let fj : M>o ^ M 6e functions satisfying IljeS fj ^ -^^"'^(I^>o) for every 
nonempty subset B C {1, . . . , A;}. Then 

j=l n=l P&V{k) BeP jeB 

Remark 3. When k = 1 Proposition [3] gives 



00 -| n fyQ 

n=l •'^ 



which agrees with [71 p. 28 (3.18)]. 

Proof. It follows from Campbell's Theorem p. 28]; cf. also [3 (3.14), (3.18) 
and Cor. 3.1]) that 

(9) e( 9l{TnMTn,)---9v{Tn^)) 

ni=nj ■^i=j 
00 



p 00 

J{^{Y.9ATn))=2-^J[ g,{x)dx 



j=l n=l j=l ' 

for all V € Z+ and aW. gi, . . . ,gv G L^(M>o). Now the left hand side of ([8]) can be 
expressed as 



E( h{Tn,)f2{Tn,)---fk{Tn, 
P£V{k) ni,...,nfeeZ+ 

where ~p is the equivalence relation on {1, . . . , A:} corresponding to P. Hence ([5]) 
follows from ^ applied for all f = 1, . . . , fc and with functions fs = rTjes /? ' B ^ P, 
in place oi gi, . . . ,gv □ 

As in Section [3] we let > 1 and fix < Vi < V2 < • • • < . We will apply 
Proposition[3]to the situation where fj = Xji ^ ^ j 1^ k, where Xj is the characteristic 
function of the interval [0, V^]. In this case the proposition states that 

k k CO 



E( 

' j=l ' ' j=l ' n=l 



(10) = 2-*''n(/ n^^(^)^^ 

p&V{k) BeP •'^ jeB 



PeV{k) BeP 
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where Jb = minj. 

Remark 4. The formula in (jlOp can also be obtained from the mixed cumulants of 
N{Vi), NiVk) (see e.g. U Prop. 6.16]). 

We are interested in comparing the formula in (|1U|) with the result of Theorem [3J 
We first need the following observation. 

Lemma 3. Let T>{k) he the set of matrices D which have nonnegative entries and 
are on the form occurring in Theorem O together with the k x k indentity matrix Ik ■ 
When D = Ik let Vi = i, 1 < i < k. Then there is a bijection g : T>(k) — ?■ V{k) with 
the property that if D ^ is an m x k matrix and g{D) = P = {-Bi, . . . , B^p} 

then #P = m and {vi, . . . ,Um} = {jBi, • • • ,3bJ\- 

Proof. For D € T>{k) of size m x k we set Bi = {j \ dij 7^ 0}, 1 < i < m, 
and define g{D) = {Bi, . . . , Bm}- It is clear that 5 is a bijection with the desired 
properties. □ 

This lemma together with Theorem [3] and the result in (jlOp immediately implies 
the following theorem. 

Theorem 4. Let Nj{L) = ^Nj{L) so that Nj{L) denotes the number of pairs of non- 
zero lattice points of L belonging to the n-hall centered at the origin having volume 
Vj, 1 < j < k. Then 

k k 

as n —7- 00. 

Corollary 1. Let k > 1 and let V = (Vi, V2, . . . , Vfc) for some fixed < Vi < 
V2 . . . < Vfe . Consider the random vectors 

N{L,V) = {Ni{L),...,Nk{L)) 

and 

N{V) = {N{Vi),...,N{Vk)). 
Then N{L,V) converges in distribution to N{V) as n ^ 00. 

Proof. Since N{Vj), I < j < k, is Poisson distributed with expectation value 
and the Poisson distribution is uniquely determined by its moments, it follows from 
[SI Thm. 3] that also the distribution of N{V) is uniquely determined by its (joint) 
moments. The corollary now follows from this, Theorem U] and a (in principle word 
by word) generalization to the present situation of [2, Thm. 30.2]. □ 

Remark 5. When k = 1 Corollary [1] shows that Ni{L) has a distribution which 
converges weakly to the Poisson distribution with mean ^Vi as n ^ 00. This was 
first proved by Rogers ([11', Thm. 3]). 

If we write A^Vi (B) instead of A^i (L) and then let Vi = t we get a stochastic process 
{Nt{L),t > 0}. It follows from Corollary [1] that all finite dimensional distributions 
coming from this process converge to the corresponding finite dimensional distribu- 
tions of the Poisson process {N(t),t > 0} as n — )• 00. Hence, by [21 Thm. 12.6 and 
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Thm. 16.7], the process {Nt{L),t > 0} converges weakly to the process {N{t),t > 0} 
as n — 7- oo. This concludes the proof of Theorem [TJ 

Remark 6. The proof of Theorem [1] can be slightly simplified by considering factorial 
moments instead of ordinary moments (cf. [5, Thm. 6.10]). However the present set- 
up will be useful also when dealing with the moments of the Epstein zeta function, 
cf. [I3j. 

Acknowledgement. The author is grateful to Andreas Strombergsson for suggesting 
the problem and for many helpful discussions and to Svante Janson for valuable 
comments. 

References 

K. Ball, A lower bound for the optimal density of lattice packings, Internat. Math. Res. Notices 
1992, no. 10, 217-221. 

P. Billingsley, Probability and measure, third edition, Wiley Series in Probability and Mathe- 
matical Statistics, John Wiley & Sons Inc., New York, 1995. 

P. Billingsley, Convergence of probability measures, second edition, Wiley Series in Probability 
and Statistics, John Wiley & Sons Inc., New York, 1999. 

L. K. Hua, Introduction to number theory. Translated from the Chinese by Peter Shiu, Springer- 
Verlag, 1982. 

S. Janson, T. Luczak, A. Rucihski, Random graphs, Wiley-Interscience Series in Discrete Math- 
ematics and Optimization, John Wiley & Sons Inc., New York, 2000. 

G. A. Kabatyanskii, V. I. Levenshtein, Bounds for packings on a sphere and in space, Problems 
of Information Transmission 14 (1978), no. 1, 1-17. 

J. F. C. Kingman, Poisson processes, Oxford Studies in Probability, vol. 3, The Clarendon 
Press, Oxford University Press, New York, 1993. 

L. C. Petersen, On the relation between the multidimensional moment problem and the one- 
dimensional moment problem. Math. Scand. 51 (1982), no. 2, 361-366. 
C. A. Rogers, Mean values over the space of lattices. Acta Math. 94 (1955), 249-287. 
C. A. Rogers, The moments of the number of points of a lattice in a bounded set, Phil. Trans. 
R. Soc. Lond. A. 248 (1955), 225-251. 

C. A. Rogers, The number of lattice points in a set, Proc. London Math. Soc. (3) 6 (1956), 
305-320. 

P. Sarnak, Values at integers of binary quadratic forms. Harmonic analysis and number theory 
(Montreal, 1996, CMS Conf. Proc, 21), Amer. Math. Soc, Providence, RI, 1997, pp. 181-203. 
W. M. Schmidt, Masstheorie in der Geometrie der Zahlen, Acta Math. 102 (1959), 159-224. 
A. Sodergren, On the value distribution and moments of the Epstein zeta function to the right 
of the critical strip, in preparation. 
[15] J. M. VanderKam, Correlations of eigenvalues on multi-dimensional flat tori. Comm. Math. 
Phys. 210 (2000), no. 1, 203-223. 

Department of Mathematics, Uppsala University, Box 480, 
SE-75106 Uppsala, Sweden 
sodergrenOmath .uu.se 



[1] 

[2] 

[3] 

[4] 

[5] 

[6] 

[7] 

[8] 

[9] 
[10] 

[11] 
[12] 

[13] 
[14] 



